Spin Waves in the Ferromagnetic Ground State 
of the Kagome Staircase System C03V2O8 
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Inelastic neutron scattering measurements were performed on single crystal C03V2O8 wherein 
magnetic cobalt ions reside on distinct spine and cross-tie sites within kagome staircase planes. 
This system displays a rich magnetic phase diagram which culminates in a ferromagnetic ground 
state below Tc ~ 6 K. We have studied the low-lying magnetic excitations in this phase within the 
kagome plane. Despite the complexity of the system at higher temperatures, linear spin- wave theory 
describes most of the quantitative detail of the inelastic neutron measurements. Our results show 
two spin-wave branches, the higher energy of which displays finite spin-wave lifetimes well below 
Tc, and negligible magnetic exchange coupling between Co moments on the spine sites. 

PACS numbers: 75.30.Ds, 75.50.Dd, 75.10.Dg 
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Magnetic materials in which the constituent magnetic 
moments reside on networks of triangles and tetrahedra 
have been of great interest due to their propensity for 
exotic ground states, a consequence of geometrical frus- 
tration While ferromagnetically-coupled moments 
on such lattices generally do not result in such ground 
states, ferromagnets, and materials which display both 
ferromagnetic (FM) and antiferromagnetic (AFM) inter- 
actions, on such lattices remain of great interest, in part 
due to the relative scarcity of well-studied examples, and 
in part due to intriguing spin ice Q and multiferroic 
phenomena 0] which characterize some of these ground 
states. 

The kagome lattice is comprised of a two-dimensional 
network of corner-sharing triangles. Several realiza- 
tions of magnetic moments on stacked kagome lattices 
with varying degrees of crystalline order have been ex- 
tensively studied. Recently studied examples include 
jarosites, such as KFe3(OH)6 (804)2 0] and herbert- 
smithite ZnCu3(OH)6Cl2 0, both of which show evi- 
dence of strong magnetic frustration. 

The stacked kagome staircase materials M3V2O8 
(M=Ni,Co) display orthorhombic crystal structures with 
space group Cmce Q. Their kagome layers are buckled 
and composed of edge-sharing M^+Og octahedra. These 
layers are separated by non-magnetic V^+04 tetrahedra. 
The buckled kagome layers are perpendicular to the or- 
thorhombic 6-axis and form what is known as a stacked 
kagome staircase structure. Figure [T] shows the projec- 
tion of the kagome staircase onto the a-c plane. The two 
inequivalent M sites are known as spines (Ml) and cross- 
ties (M2). The superexchange interaction between spine 
and cross-tie sites and between two adjacent spine sites 
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FIG. 1: [color online] (a) A schematic diagram of the kagome 
staircase structure as reduced to 2D in the a-c plane. The 
cobalt ions are represented by open and solid circles for spine 
(Ml) and cross-tie sites (M2), respectively. Chains of spine 
sites running parallel to the a-direction are alternatively above 
(0) and below the plane (®). (b) The basis used in the linear 
spin-wave theory calculation. 



are denoted by Jgc and Jss, respectively. 

One member of this family, Ni3V208 (NVO), under- 
goes a series of phase transitions on lowering tempera- 
ture very interesting characteristic of 
this compound is that it exhibits simultaneous ferroelec- 
tric and incommensurate AFM order, that is, multifer- 
roic behaviour, in one of its ordered phases. In isostruc- 
tural C03V2O8 (CVO), the S — 1 magnetic moments at 
the Ni^+ site are replaced with S — 3/2 Co^+ moments. 
CVO also displays a rich low temperature phase diagram, 
which has been studied using polarized and unpolarized 
neutron diffraction, dielectric measurements [l2| , magne- 
tization and specific heat measurements 13|. There is a 



series of four AFM ordered phases below Tm = 11-3 K 
which can be characterized by incommensurate or com- 
mensurate ordering wavevectors (OrO). In contrast to 
NVO, the ultimate ground state in CVO is ferromagnetic 
and the Curie temperature is Tc ~ 6 K. Earlier powder 
neutron diffraction measurements on CVO showed 
ordered magnetic moments of 2.73(2) and 1.54(6) /is on 
the spine and cross-tie sites, respectively, at 3.1 K. All 
moments are aligned along the a-axis direction. 

While much work has been performed on the phase di- 
agrams of NVO and CVO, little is known about the exci- 
tations and, correspondingly, the underlying microscopic 
spin Hamiltonian for these topical magnets. The ferro- 
magnetic state in CVO is an excellent venue for such a 
study, as the ground state itself is very simple and there- 
fore the excitations out of the ground state should be 
amenable to modeling. In this Letter we report inelastic 
neutron scattering (INS) measurements of the spin- wave 
excitations within the kagome staircase plane in the FM 
ground state of single crystal CVO. These measurements 
are compared with linear spin wave theory which shows 
a surprising sublattice dependence to the exchange inter- 
actions. 

A large (5 g) and high-quality single crystal of CVO 
was grown using an optical floating-zone image furnace 
[l3 |. Thermal INS measurements were performed at the 
Canadian Neutron Beam Centre at the Chalk River Lab- 
oratories using the C5 triple-axis spectrometer. A py- 
rolytic graphite (PC) vertically-focusing monochroma- 
tor and flat analyzer were used. Measurements were 
performed with a fixed final neutron energy of Ef = 
13.7 meV and a PC filter in the scattered beam. The 
collimation after the monochromator was 29'-34'-72' re- 
sulting in an energy resolution of 0.9 meV FWHM. 

The crystal was oriented with the {hO£) kagome 
staircase plane coincident with the scattering plane. 
Constant-Q energy scans were performed along the high 
symmetry (/lOO) and (OOi) directions in this plane. Fig- 
ures 2 (a) and (b) show representative constant-Q scans 
at T = 3 K and Q=(1.8, 0, 0) and (0, 0, 1.8), respec- 
tively. The (002) nuclear Bragg peak is very weak, and is 
coincident with a strong FM Bragg peak below Tc- The 
inset of Fig. [2lja) shows the temperature dependence of 
this Bragg reflection on independent warming and cool- 
ing runs. An abrupt falloff in intensity near Tc ~ 6 K 
and accompanying hysteresis indicate the strongly dis- 
continuous nature of this phase transition. 

The overall features of the two spectra in Fig. 2 are 
quantitatively similar at T = 3 K. Two spin-wave exci- 
tations, identified on the basis of their temperature and 
Q-dependencies, are observed and have been modelled 
using resolution-convoluted damped harmonic oscillator 
(DHO) lineshapes IJ]. The resulting fits are shown as 
the dashed lines in Fig. 2, and this analysis allows us to 
conclude that the higher-energy spin wave, near 5.7 meV 
in both cases, has a substantial intrinsic energy width 
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FIG. 2: Representative constant-Q INS spectra with (a) Q= 
(1.8,0,0) and (b) Q= (0,0, 1.8) at T = 3 K (FM phase). The 
broken line shows the resolution-corrected fits to the data, as 
described in the text. The solid horizontal bars indicate the 
instrumental energy resolution. The inset shows heating and 
cooling scans of the (mainly) magnetic elastic (002) Bragg 
scattering, characteristic of the FM ground state in CVO. 



of r=0.70(8) meV at Q=(1.8,0,0) and 0.90(8) meV at 
Q=(0, 0,1.8). 

A series of constant-Q scans for {hOO) and (00£) direc- 
tions were collected at T = 3 K and are presented as a 
color contour map in Figs. 3 (a) and (c). Dispersive fea- 
tures corresponding to two bands of spin waves are seen 
in both data sets. The top of the upper spin- wave band at 
AE ^5.7 meV corresponds to excitations reported ear- 
lier using a time-of-flight technique Q . These constant-Q 
scans were fit to resolution-convoluted DHO lineshapes, 
which gave intrinsic energy widths for the higher-energy 
spin-wave mode at all wavevectors ranging from F = 0.6 
to 1.1 meV, while the lower-energy spin waves were res- 
olution limited at all wavevectors. This indicates a finite 
lifetime for the higher-energy spin waves even at temper- 
atures - Tc/2. 

The spin-wave spectrum evolves rapidly with increas- 
ing temperature. Figured (a) shows Q = (3.4, 0, 0) scans 
at T = 3 K (FM phase) and T 20 K [well into the para- 
magnetic (PM) phase] . In the FM phase there is a promi- 
nent spin-wave peak at AE = 2.0 meV. A higher energy 
spin-wave peak is also present in this scan but with a 
much lower peak intensity and an intrinsic energy width 
of r = 0.9 meV. At 20 K the well-defined low-energy spin 
wave has disappeared and only broad low-AiJ scattering 
remains. The inset of Fig. [4] (a) shows the temperature 
dependence of spin wave peak intensity at AE = 2.0 meV 
in the neighborhood of Tc- The same rapid falloff as was 
seen in the magnetic Bragg scattering at (002) is seen 
in the inelastic intensity, as well as the same hysteresis. 
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FIG. 3: Color contour maps of INS at T = 3 K [(a) and 
(c)] and corresponding linear spin-wave theory [(b) and (d)] 
as described in the text. The broken lines show the dispersion 
relations resulting from this spin-wave theory analysis. 
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FIG. 4: [color online] (a) INS scans at Q = (3.4, 0, 0) for 
T = 3 K and T = 20 K. The inset shows the temperature 
dependence of the scattering at the inelastic peak position. 
The lines are guides to the eye. (b) Scans in the PM phase 
for Q = (1.8,0,0) and Q = (0,0,1.8). The inset (same leg- 
end as main figure) shows the temperature dependence of the 
inelastic scattering at the spin-wave peak (for T — 3 K, see 
Fig.E]) at AE = 1.86 meV for the different Q directions. The 
vertical dotted lines in the inset indicate Tc and Tjv. 



indicating that these spin waves are strongly coupled to 
the ferromagnetic order parameter. 

Figure 2] (b) shows the same constant-Q scans as in 
Fig. m but now at 20 K rather than 3 K. The well- 
defined spin- wave modes are no longer present, and the 
low-energy inelastic scattering is significantly greater at 
Q=(0, 0,1.8) as compared with (1.8,0,0) at T = 20 K. 
The temperature dependence of this scattering is con- 
trasted in the inset of Fig. H (b). The Q=(1.8,0,0) in- 
elastic scattering drops quickly at Tc ~ 6 K while that at 
Q= (0, 0, 1.8) gradually increases with temperature. We 
attribute this difference to prominent longitudinal easy- 
axis spin fluctuations along the a-axis at high temper- 
atures. Given that the neutron scattering cross-section 
IJj is sensitive to magnetic fluctuations perpendicular 
to Q, longitudinal fluctuations would be seen in the 
Q=(0, 0,1.8) spectrum rather than the (1.8,0,0) spec- 
trum. Note that the highest temperature transition, to 
the paramagnetic state at Tn = 11.3 K, is evident as a 
change in slope of the temperature dependence shown 
in the inset of Fig. [4] (b), and both Tc ~ 6 K and 
T/v = 11.3 K are indicated by dashed vertical lines in 
this inset. 

We have carried out a linear spin-wave theory analysis 
of the magnetic excitations observed in Figs. 3 (a) and 
(c) to understand as much of the relevant microscopic 
spin Hamiltonian as possible. The full Hamiltonian is 
potentially complicated if account is taken of the two in- 



equivalent magnetic sites and the 3D kagome staircase 
structure. We employed a 2D model in which the mag- 
netic ions in a layer are projected onto the average plane 
of the layer (Fig.[T|) and only near-neighbor exchange and 
single-ion anisotropy are included. 

We change from the conventional centered-rectangular 
unit cell to a primitive rhombohedral unit cell defined 
by lattice vectors Ri = (§,§) and R2 = (§,— f) and 
basis vectors <5i — Ri/2, ^2 — R2/2, and 63 = 0. If 
R describes the set of lattice points and Sr is the spin 
operator at a location r we can write 



Ti-sc = ^Jsc'^^ ^ Sr-Sr+5, (1) 

R S=±Si,±d2 

Ti-ss = — ^ss ^ [Sr+6i -Sr-^^ + Sr_5j -Sr+s^] . (2) 
R 



Ti.sc and Tigs are the exchange interactions between spine 
and cross-tie spins with couplings Jsc and Jss- The fact 
that the spine and cross-tie spins are found to be fer- 
romagnetically aligned implies that Jsc is positive. We 
choose the spin z-axis parallel to the crystallographic 
a-axis, consistent with the ordered moment direction. 
Magnetization measurements have established that the 
crystallographic 6-axis (spin x-axis) is a harder axis than 
the c-axis (spin y-axis) [13|. The appropriate single-ion 
anisotropy Hamiltonian, Tia, distinguishes the three or- 
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thogonal directions and the inequivalent sites 

^'^ = E E E ^fis^+sf- (3) 

R i=l,2,3 a=x,y,z 

The sequence from hard to easy axis is estabhshed by the 
condition > A\ > Af. We use reduced anisotropics: 



A, = (Af +Ay- 2Al)/2 and = (Af - A\)/2. Since 
the i = 1 and i = 2 sites are equivalent there are 4 
independent parameters A^, Ac, e^, and ec- 

In hnear spin- wave theory ^15^ the total Hamiltonian 
Ti may be written in term of spin- wave operators Ci(k) 
and c|(k) 



W = ^ ;i.,(k)4(k)c,(k) + 1 ^ 5,,(k) (4(k)c](-k) + c,(-k)c,(k) 



(4) 



f^SC H" J Si 



A,, 



/i(k) = 2S'( -Jsscosk- (<5i + ^2) 
— Jsc cos k • <5i 



— Jss cos k • (^1 + (52) — Jsc cos k-^i 



J sc ^" Jss ^ As 
— Jsc cos k • ^2 



— Jsc cos k • ^2 

2 Cf \ 



, 

5(k) = 25 I e, 




(5) 



where i,j — 1, 2, 3. Although the moments on the differ- 
ent sites are unequal we make the simplifying assumption 
that S = 3/2 for all sites. 

The unpolarized magnetic neutron scattering cross- 
section and the spin-spin correlation functions it contains 
can be related to one spin-wave Green's functions. The 
Green's functions can be calculated by inverting a matrix 
involving the two matrices, /i(k) and (/(k) defined above, 
following Coll and Harris [l6|. The spin- wave energies 
are determined by solving for the zeros of a determinant 
involving matrix elements of the inverse Green's func- 
tion. The resulting scattering function S{Q, AE) for Q 
parallel to (00^) is proportional to 

id 

where / is the identity matrix, k = Q and z = AE — iT 
is the complex energy, and i and j are summed over the 
indices of the 3x3 matrix in square brackets. The ex- 
pression for Q parallel to (hOQ) has a similar form. The 
resulting scattering functions are multiplied by the mag- 
netic form factor, the Bose factor, and a single intensity 
scale factor, and are plotted in Fig. [3] (b) and (d). 

Best agreement between the experimental data and the 
spin-wave theory calculation in Fig. [3] was obtained for 
magnetic coupling predominantly between the spine and 
cross-tie Co ions with Jsc = 1.25 ± 0.08 meV, while the 
spine-spine coupling Jss vanishes. The best fit spin-wave 
uniaxial anisotropy parameters are As = 1.5 ± 0.1 meV, 
Ac = 2.13 ± 0.2 meV, = 0.6 ± 0.3 meV, and ec = 
1.2 ±0.3 meV. 

Figure [3] shows that the spin- wave theory gives a very 
good description of the dispersion of the two modes 
(dashed lines) and accounts for the observed trend of 
the spin waves to trade intensity as a function of Q. 
This description is not perfect, however. The calculated 



dispersion of the lower spin-wave band is low compared 
with experiment near (200) and (002) where the inten- 
sity is very weak. The calculation is not convolved with 
the instrumental resolution; instead the energy width is 
manually set in both high and low-energy bands to corre- 
spond to the measured width. The broad (in energy) neu- 
tron groups corresponding to the upper spin- wave bands 
are most evident near (200) and (002). The lower en- 
ergy spin-wave band becomes much more intense near 
the zone centers of (400) and (004). Steep excitation 
branches near (400) and (004) with comparitively weak 
intensity in the experimental data [Figs. 3 (a) and (c)] 
are identified as acoustic phonons with a speed of sound 
of 1050 ± 100 m/s, in both directions. 

To conclude, our INS study of the spin-wave excita- 
tions in the ferromagnetic ground state of CVO within 
its kagome staircase plane reveal two separate spin-wave 
bands between 1.6 and 5.7 meV. The upper spin-wave 
band is damped with finite energy widths F in the range 
of 0.6 to 1.1 meV. These spin-wave excitations can be 
accurately described by a simple model Hamiltonian and 
linear spin-wave theory. The model gives a magnetic cou- 
pling that is predominantly between the spine and cross- 
tie sites of the kagome staircase. 
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